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Abstract 

We present an improved measurement of the standard model CP violation parameter sin2<^i 
(also known as sin 2(3) based on a sample of 85 x 10 6 BB pairs collected at the T(45) resonance 
with the Belle detector at the KEKB asymmetric-energy e + e~ collider. One neutral B meson 
is reconstructed in a J/ipK^, ip(2S)Kg, XciK^ VcKs^ J/i>K*°, or J/^Kl CP-eigenstate decay 
channel and the flavor of accompanying B meson is identified from its decay products. From the 
asymmetry in the distribution of the time intervals between the two B meson decay points, we 
obtain sin2</>i = 0.719 ± 0.074(stat) ± 0.035(syst). We also report measurements of CP violation 
parameters for the related B° — > J/t/jir decay mode and the penguin-dominated processes B° — > 
n'K° s , (pKg and K + K-R° s . 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw 
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In the Standard Model (SM), CP violation arises from an irreducible complex phase in 
the weak interaction quark- mixing matrix (CKM matrix) In particular, the SM predicts 
a CP- violating asymmetry in the time-dependent rates for B° and B° decays to a common 
CP eigenstate, fcp, with negligible corrections from strong interactions 0: 

A f .\ _ ^(-B — »• f C p) — V(B — >• fcp) t . , . , A » m 
Mt) = r(BO->fcp) + T(B°^f C p) = ~^ 2 ^ Am ^ C 1 ) 

where T(B°,B° — > /cp) is the decay rate for a 5° or 5° to fcp dominated by a 6 — > ccs 
transition at a proper time t after production, £/ is the CP eigenvalue of fcp, Arrid is the 
mass difference between the two B° mass eigenstates, and <\>\ is one of the three interior 
angles of the CKM unitarity triangle, defined as <pi = n — arg(— V^Vu/ ~ V^Vcd)- Non-zero 
values for sin20i were reported by the Belle and BaBar groups^ 01. 

Belle's published measurement of sin20i is based on a 29.1 fb _1 data sample containing 
31.3 x 10 6 BB pairs produced at the T(4S) resonance. In this paper, we report an improved 
measurement that uses 85 x 10 6 BB pairs (78 fb _1 ). The data were collected with the Belle 
detector || at the KEKB asymmetric collider ||, which collides 8.0 GeV e~ on 3.5 GeV e + 
at a small (±11 mrad) crossing angle. We use events where one of the B mesons decays to 
fcp at time top, and the other decays to a self-tagging state, /tag, i-e., a final state that 
distinguishes B° and B°, at time tt ag - The CP violation manifests itself as an asymmetry 
A(At), where At is the proper time interval between the two decays: At = tcp — ttag- 
At KEKB, the T(4S f ) resonance is produced with a boost of fij = 0.425 nearly along the 
electron beam direction (z direction), and At can be determined as At ~ Az/{f3^f)c, where 
Az is the z distance between the fcp and / tag decay vertices, Az = zqp — z tag . The Az 
average value is approximately 200 /im. 

The Belle detector [BJ is a large-solid-angle spectrometer that consists of a silicon vertex 
detector (SVD), a central drift chamber (CDC), an array of aerogel threshold Cerenkov coun- 
ters (ACC), time-of-flight scintillation counters (TOF), and an electromagnetic calorimeter 
comprised of CsI(Tl) crystals (ECL) located inside a superconducting solenoid coil that pro- 
vides a 1.5 T magnetic field. An iron flux-return located outside of the coil is instrumented 
to detect K\ mesons and to identify muons (KLM). 

We reconstruct B° decays to the following CP eigenstates 0: J/ipKg, ip(2S)Kg, XciKg, 
rj c Kg for £f = —1 and J/ipK^ for = +1. We also use B° — > J/tpK*° decays where 
K*° — > Kg7t°. Here the final state is a mixture of even and odd CP, depending on the 
relative orbital angular momentum of the J/ip and K*°. We find that the final state is 
primarily £y = +1; the £f = —1 fraction is 0.19±0.02(stat)±0.03(syst)||. For reconstructed 
B — ► fcp candidates other than J/i/jK®, we identify B decays using the energy difference 



AE = E^ 1S — i?b™ s m and the beam-energy constrained mass M bc = yj (E^ m ) 2 — (pg 118 ) 2 , 
where E^™ m is the beam energy in the center-of-mass system (cms), and E c ™ and p c ^ s are 
the cms energy and momentum of the reconstructed B candidate, respectively. Figure |l| 
(left) shows the M^ c distributions for all B° candidates except for B° — > J/ipK^ that have 
AE values in the signal region. Table | lists the numbers of observed candidates (iV rec ). 

Candidate B° — ► J/ipK^ decays are selected by requiring ECL and/or KLM hit patterns 
that are consistent with the presence of a shower induced by a neutral hadron. The centroid 
of the shower is required to be in a 45° cone centered on the K\ direction that is inferred from 
two-body decay kinematics and the measured four-momentum of the J/ip. Figure |I| (right) 
shows the p c ™ distribution, calculated with the B° — ► J/i/jK® two-body decay hypothesis. 
The histograms are the results of a fit to the signal and background distributions. There 
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FIG. 1: The beam-energy constrained mass distribution for all decay modes combined other than 
J/tpK^ (left). The p c ^ s distribution for B° — > J/ipK^ candidates with the results of the fit (right). 



TABLE I: The numbers of reconstructed B — > fcp candidates before flavor tagging and vertex 
reconstruction (N iec ), the numbers of events used for the sin2(/>i determination (iV ev ), and the 
estimated signal purity in the signal region for each fcp mode. 



Mode 


if 


iVrec 


N cv 


Purity 


J/^(£ + £-)K° s (tt + tt-) 

J/4>(l + l-)K° s {iT Tr ) 

^{2S)(£ + e-)K° s (7T + 7T-) 

ifj(2S)(J/iJ7r + 7r-)K s (7r+ir-) 

Xci(J/^)K s (it+tt-) 


-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 


1285 
188 
91 
112 
77 
72 
49 
21 


1116 
162 
76 
96 
67 
63 
44 
15 


0.98 
0.82 
0.96 
0.91 
0.96 
0.65 
0.72 
0.94 


All with £f = -1 




1895 


1639 


0.94 


J/^(e + r)K*°(K s 7r°) 


-1(19%)/+1(81%) 


101 


89 


0.92 


j/^(e + £-)K° L 


+1 


1330 


1230 


0.63 


All 




3326 


2958 


0.81 



are 1330 entries in total in the 0.20 < p§° s < 0.45 GeV/c signal region; the fit indicates a 
signal purity of 63%. The reconstruction and selection criteria for all of fcp channels used 
in the measurement are described in more detail elsewhere 0. 

Leptons, charged pions, kaons, and A baryons that are not associated with a reconstructed 
CP eigenstate decay are used to identify the 6-flavor of the accompanying B meson: high 
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momentum leptons from b — > c£~T7; lower momentum leptons from c — > s£ + i>; charged 
kaons and A baryons from b — > c — > s; fast pions from B° — > (n + , p + , af , etc.); and 

slow pions from D*~ — > D°tt~. Based on the measured properties of these tracks, two 
parameters, q and r, are assigned to an event. The first, q, has the discrete values q = ±1 
that is +1 (—1) when B ta . g is likely to be a B° (B°), and the parameter r is an event-by-event 
Monte-Carlo-determined flavor-tagging dilution factor that ranges from r = for no flavor 
discrimination to r = 1 for an unambiguous flavor assignment. It is used only to sort data 
into six intervals of r, according to flavor purity; the wrong-tag probabilities, w% (I = 1,6), 
that are used in the final fit are determined directly from data. Samples of B° decays to 
exclusively reconstructed self-tagged channels are utilized to obtain w\ using time-dependent 
B°-B° mixing oscillation: (Aof — A^sf)/(Aof + Asf) = (1 — 2w{) cos(AmdAt), where Nop 
and iVsF are the numbers of opposite and same flavor events. The event fractions and wrong 
tag fractions are summarized in Table |T[ The total effective tagging efficiency is determined 
to be X)f=i //(I — 2w/) 2 = 0.288 ± 0.006, where fi is the event fraction for each r interval. 

TABLE II: The event fractions (ej) and wrong tag fractions (wi) fo each r interval. The errors 
include both statistical and systematic uncertainties. The event fractions are obtained from the 
J/ipKg simulation. 



I r ei wi 

1 0.000 - 0.250 0.399 0.458 ± 0.006 

2 0.250 - 0.500 0.146 0.336 ± 0.009 

3 0.500 - 0.625 0.104 0.229 tomi 

4 0.625 - 0.750 0.122 0.159 ±0.009 

5 0.750 - 0.875 0.094 0.111 ±0.009 

6 0.875 - 1.000 0.137 0.020 ±2*2SI 



The vertex position for the fcp decay is reconstructed using leptons from J/xjj decays 
or kaons and pions from r\ c and that for f tag is obtained with well reconstructed tracks 
that are not assigned to fcp- Tracks that are consistent with coming from a Kg — > 7t + tt~ 
decay are not used. Each vertex position is required to be consistent with a run-by-run- 
determined interaction region profile that is smeared in the r-0 plane by the B meson decay 
length. With these requirements, we are able to determine a vertex even with a single track; 
the fraction of single-track vertices is about 10% for zqp and 22% for z tag - The proper-time 
interval resolution function, i? sig (At), is formed by convolving four components: the detector 
resolutions for zqp and z t ag, the shift in the z taug vertex position due to secondary tracks 
originating from charmed particle decays, and smearing due to the kinematic approximation 
used to convert Az to At. A small component of broad outliers in the Az distribution, 
caused by mis-reconstruction, is represented by a Gaussian function. We determine twelve 
resolution parameters from the data from fits to the neutral and charged B meson lifetimes 
and obtain an average At resolution of ~ 1.43 ps (rms). The width of the outlier component 
is determined to be (42l|) ps; the fractional areas are (2 ± 1) x 10~ 4 and (2.7 ± 0.2) x 10~ 2 
for the multiple- and single-track cases, respectively. 

After flavor tagging and vertexing, we find 1465 events with q = ±1 flavor tags and 
1493 events with q — — 1. Table | lists the numbers of candidates used for the sin20i 
determination (N ev ) and the estimated signal purity in the signal region for each fcp mode. 
Figure |2] shows the observed At distributions for the q£f = ±1 (solid points) and q£f = — 1 
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(open points) event samples. The asymmetry between the two distributions demonstrates 
the violation of CP symmetry. We determine sin 20i from an unbinned maximum-likelihood 




At (ps) 



FIG. 2: At distributions for the events with q£f = +1 (solid points) and q£f = — 1 (open points). 
The results of the global fit with sin 2</>i = 0.719 are shown as solid and dashed curves, respectively. 

fit to the observed At distributions. The probability density function (pdf) expected for the 
signal distribution is given by 

e -|Ai|/r s0 

V^(At,q,wi,^ f ) = — [1 -q£ f (l -2wi)sm2<j> 1 sin(Am d At)], (2) 

At b o 



where we fix the B° lifetime (rgo) and mass difference at their world average values [10 



Each pdf is convolved with the appropriate R S i g (At) to determine the likelihood value for 
each event SIS db function of sin 2<pi\ 

Pi = (1 - /oi) / [/sigTWAf, q, w h Z f )R sig (At - At') 

+ (1 - f sig )V hkg (At')R bkg (At - At')]dAt' + / iP ol (At), (3) 

where / s i g is the signal probability calculated as a function of p ™ for J/i/jK® and of AE 
and Mbc for other modes. 7\kg(At) is the pdf for combinatorial background events, which 
is modeled as a sum of exponential and prompt components. It is convolved with a sum 
of two Gaussians, -Rbkg, which is regarded as a resolution function for the background. To 
account for a small number of events that give large At in both the signal and background, 
we introduce the pdf, P Q \, and the fractional area, / Q i, of the outlier component. The only 
free parameter in the final fit is sin20i, which is determined by maximizing the likelihood 
function L = fli Pi, where the product is over all events. The result of the fit is 

sin 20x = 0.719 ± 0.074(stat) ±0.035(syst). 

The systematic error is dominated by uncertainties in the vertex reconstruction (0.022). 
Other significant contributions come from uncertainties in: Wi (0.015); the resolution func- 
tion parameters (0.014); a possible bias in the sin20i fit (0.011); and the J/ipK® background 
fraction (0.010). The errors introduced by uncertainties in Amj and r^o are less than 0.010. 
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FIG. 3: (a) The raw asymmetry for all modes combined. The asymmetry for J/t/jK^ and J/ipK*° 
is inverted to account for the opposite CP eigenvalue. The corresponding plots for (b) (cc)Kg, 
(c) J/ipK®, and (d) non-CP control samples are also shown. The curves are the results of the 
unbinned maximum likelihood fit applied separately to the individual data samples. 
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TABLE III: The numbers of candidate events (iV ev ) and values of sin2(/)i for various subsamples 
(statistical errors only). 



sample 


A cv 


sin 2(f>i 


J /tpKg{7T^7T ) 


111b 


0.73 ± 0.1L) 


(rc\K® PYPPnt J 1 lil TC^ ( TC~^~ 7T ^ 


593 




J/^Kl 


1230 


0.78 ±0.17 


J/i>K*°{K° s 7r ) 


89 


0.04 ± 0.63 


/ tag =B°(q = +1) 


1465 


0.65 ±0.12 


/ tag = 5° (g = -l) 


1493 


0.77 ±0.09 


< r < 0.5 


1600 


1.26 ±0.36 


0.5 < r < 0.75 


658 


0.62 ±0.15 


0.75 < r < 1 


700 


0.72 ± 0.09 


All 


2958 


0.72 ±0.07 



A number of checks on the measurement are performed. Table [TTT| lists the results obtained 
by applying the same analysis to various subsamples. All values are statistically consistent 
with each other. Figure ^(a), (b), and (c) show the raw asymmetries and the fit results for 
all modes combined, (cc)Kg, and J/vpK®, respectively. A fit to the non-CP eigenstate self- 
tagged modes B° — > D^~tt + , D*~p + , J/ipK*°(K + 7i~), and D*~£ + v, where no asymmetry 
is expected, yields 0.004 ± 0.015(stat). Figure |||(d) shows the raw asymmetry for these 
non-CP control samples. 

The signal pdf for a neutral B meson decaying into a CP eigenstate (Eq. (H)) can be 
expressed in a more general form as 



V sig (At,q,w h £f) 



-[Atl/TV 



4r B o 



2w, 



2\\\(-£ f )acp 



|A| 2 ±1 



sin(AmrfAt) 



+■ 


A 


2 - 1 




A 


2 ±1 



- cos(AmdAt) | 



(4) 



where A is a complex parameter that depends on both B°-B° mixing and on the amplitudes 
for B° and B° decay to a CP eigenstate. The parameter acp in the coefficient of sm(Am^At) 
is given by acp = — £//mA/|A| and is equal to sin20i in the SM. The presence of the cosine 
term (|A| ^ 1) would indicate direct CP violation; the value for sin20x reported above 
is determined with the assumption |A| = 1, as expected in the SM. In order to test this 
assumption, we also performed a fit using the above expression with acp and |A| as free 
parameters, keeping everything else the same. We obtain 

|A| = 0.950 ± 0.049(stat) ± 0.026(syst) 



and acp = 0.720 ± 0.074(stat) for all CP modes combined, where the sources of the sys- 
tematic error for |A| are the same as those for sin 20!. This result confirms the assumption 
used in our analysis. 

Finally, we also report on the time- dependent CP asymmetries in the B° — » rj'Kc. 
<f)Kg and B° — > K + K~Kg decays that are dominated by the b —>■ s penguin diagrams, 
and in the B° — > J/ipir decay governed by the Cabibbo-suppressed b —>■ ccd transition. 
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All measurements are based on 85 x 10 6 BB pairs. The reconstruction methods for these 
modes are described elsewhere 
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The signal yields are summarized in Table |V[ The 
flavor-tagging, vertexing and fitting techniques are the same as those described above. The 
decay rate has a time dependence given by 



V sig {At,q,wi) 



-|At|/r s0 

4r B o 



(l + q(l- 2wi)[S sm(Am d At) + .4cos(Am d At)]}. (5) 



Table [TV] summarizes the results of the fits to the At distributions. In the table, we show 
the values of — denoted by "sin 20]", since it should be equal to sin 20! in the limit that 
new physics in the penguin loop does not contribute additional CP- violating phases. 



TABLE IV: Results of the CP asymmetry parameter measurements for B° — > rj'Kg, B° — > 4>K^, 
B° — ► K + K~ Kg, and B — > J/i/jtt° decays. "sin20i" is defined as — £/<S and should be equal to 
sin 20i in the limit of no contributions from new physics in the penguin loop. For "sin20i" and 
A, the first errors are statistical and the second errors are systematic. The third errors for the 
K + K~ K s mode arise from the uncertainty in the fraction of the CP-odd component [12|. 



Mode ^ iV cv Purity (%) "sin20i" A 

7]'K S -L 299 49 +0.76 ± 036™! +0.26 ± 0.22 ± 0.03 

4>K S -1 53 67 -0.73 + 0.64 + 0.18 -0.56 + 0.41 + 0.12 

K + K~Kg —1(3%)/ + 1(97%) 191 49 +0.52 + 0.46 + 0.11^-0.42 + 0.36 + 0.09^ 
J/4>tt° +1 57 86 +0.93 + 0.49 + 0.08 -0.25 + 0.39 + 0.06 
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